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ABSTRACT 


Optimum  linear  and  quadratic  discrimination  filtering  tech- 
niques are  developed  for  discriminating  between  short  period 
seismic  records  originating  from  earthquakes  and  explosion.. 
Linear  and  quadratic  detection  riltering  and  matched  filtering 
are  compared  with  the  classical  spectral  ratio  and  complexity 
measure  using  a learning  population  of  LASA  array  beair^  of  23 
earthquakes  and  15  explosions  and  a test  population  with  17 
earthquakes  and  11  explosions.  The  linear  detection  filter 
misc lass i f ies  one  event  in  the  test  set  whereas  all  othe?  tech- 
niques misclassify  between  three  and  five  events. 

It  is  shown  that  for  the  spectral  ratio  discriminant  the  dis 
criminatory  power  lies  in  the  ratio  of  .4-. 8 Hz  energy  to  1.0  Hz 
energy,  and  that  the  higher  frequency  energy  has  no  additional 
discriminatory  power.  We  find  that  the  explosions  which  fail  to 
discriminate  were  probably  cratering  experiments.  In  either 
this  case,  or  if  the  explosion  is  deep,  pP  will  not  cancel  P at 
low  frequencies.  Thus,  we  propose  that  pP-P  cancellation  is 
the  basic  physical  explanation  for  the  success  of  short  period 
discrimination . 
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INTRODUCTION 


The  use  of  short  period  datu  in  discriminating  between  seis- 
mic records  originating  from  earthquakes  and  explosions  can  be 
important,  especially  when  a reliable  estimate  of  the  suriace 
wave  magnitude  cannot  be  obtained  from  the  corresponding  long 
period  records.  In  the  case  where  both  long-period  surface  wave 
and  short-period  body  wave  magnitudes  are  measured,  one  can  achieve 
fairly  good  discr iminab i 1 ity  on  the  basis  of  these  two  measurements 
alone  (sec,  for  example,  Lacoss,  1909;  Capon  ct  al.,  1909;  Shumway 
and  Blandford,  1970;  Shumway  and  Blandford,  1972).  However,  it 
is  reasonable  to  suppose  that  the  short-period  waveforms  contain 
more  discrimination  information  than  can  be  inferred  solely  from 
the  magnitude  of  the  body  wave.  This  speculation  is  coupled  with 
the  fact  that  the  body  wave  magnitude  is  useless  without  the 
corresponding  surface  wave  magnitude  whicn  may  not  be  observable. 
The  above  objections  justify  a more  detailed  investigation  of 
the  discriminatory  capabilities  of  the  short-period  waveforms. 

The  two  most  prominently  mentioned  spectral  discriminants 
arc  ter  ted  the  spectral  ratio  and  complexity.  The  spectral  ratio 
operates  on  the  hypothesis  that  the  ratio  of  low  frequency  energy 
to  high  frequency  energy  for  the  source  spectrum  of  earthquakes 
exceeds  that  for  explosions.  (Sec,  for  example,  Booker  and 
M i tronovas , 19b4.)  Blandford,  1974  , however,  has  developed  a 
theory  suggesting  that  this  hypothesis  docs  not  always  hold.  In 
a technical  report  Kelly,  1908,  and  two  reports  of  Lacoss,  1909, 
showed  that  the  spectral  ratio  could  produce  a reasonable  separa- 
tion when  applied  to  suites  of  earthquakes  and  explosions.  A 
second  discriminant  is  suggested  if  one  assumes  that  explosions 
are  basically  impulsive  while  earthquakes  generally  have  a more 
complicated  waveform,  A measure  of  complexity  consisting  of  the 
total  energy  in  the  signal  coda  has  been  applied  in  Carpenter, 

1904,  and  Kelly,  1968,  with  some  success. 
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An  appreciation  lor  the  possible  discriminatory  capabilities 
ot  the  two  discriminants  may  be  gained  by  inspecting  figure  1 which 
summarizes  tne  data  from  the  40  presumed  earthquakes  and  2b  pre- 
sumed explosions  given  in  Figures  A1  and  A2  in  Appendix  I.  The 
traces  are  LASA  beams  for  events  located  between  40  and  b0°.\ 
latitude  and  between  0 and  90  ° L longitude  plus  other  selected 
"difficult”  events  (see  Appendix  1).  The  original  data  traces 
were  gain  adjusted  to  a unit  amplitude,  aligned  on  first  motion 
and  averaged  to  produce  the  earthquake  and  explosion  signals  shown 
in  the  left  halt  ol  Figure  1.  lie  sec  that  the  earthquake  signal 
tends  to  have  more  energy  in  the  coda  suggesting  that  complexity 
might  be  a useful  discriminant.  The  average  earthquake  and  explo- 
sion noise  spectra  were  computed  by  subtracting  from  each  trace 
its  mean  signal  (the  mean  explosion  or  mean  earthquake)  and  then 
computing  tne  average  spectrum  of  the  residuals.  These  "error 
spectra  tend  to  display  greater  power  at  low  frequencies  for 
the  earthquakes  indicating  that  a spectral  ratio  could  be  used  as 
a discriminant.  Averages  from  a randomly  drawn  subset  of  the 
population,  shown  on  the  right  half  ot  Figure  1,  suggests  that 
the  observations  made  above  hold  fairly  consistentlv  for  both 
earthquakes  and  explosions. 

The  earthquake  and  explosion  signals  in  Figure  1 also  suggest 
that  more  information  may  be  available  than  would  be  provided  by 
the  previously  discussed  spectral  and  complexity  measures.  For 
example,  the  rather  regular  pattern  in  the  mean  earthquake  coda 
is  essentially  the  averaging  of  the  later  phases  which  can  be 
observed  on  the  original  earthquake  traces  in  Figures  A1  and  A2. 

1 he  consistency  displayed  by  the  mean  codas  of  the  full  and  subset 
populations  indicates  that  the  pattern  might  be  a general  charac- 
teristic of  earthquakes  from  this  region.  On  the  other  hand,  they 
may  be  etfccts  associated  with  the  geographical  source  regions, 
and  an  explosion  set  off  at  the  earthquake  site  might  show  the 
same  phasing.  Arguing  heur istical ly  from  Figure  1,  there  are 
reasons  for  believing  that  earthquakes  and  explosions  are  com- 
posed of  (a)  fixed  mean  waveforms  which  are  different  and 


(b)  stochastically  varying  noises  with  different  spectra.  If 
the  residual  noise  processes  are  Gaussian  and  stationary,  the 
problem  may  be  formulated  in  classical  discrimination  detection 
terms.  That  is,  denoting  the  received  waveform,  signal  and 
noise  (not  earth  noise)  re-'1  actively  by  x,  s and  n,  we  have  the 
two  hypotheses  Hj  (Earthquake)  and  H-,  (Explr  ion)  specifying 

H . : x = s , + n 

1 J. 

h2:  x = s 2 + n (1.1) 

where  n has  zero  means  and  autocorrelation  matrices  and  R2 
respectively  under  and  W^.  In  this  case  the  optimum  detector 
is  the  sum  of  the  outputs  of  a quadratic  filter  and  a linear 
filter,  where  optimality  is  defined  as  maximizing  the  explosion 
detection  probability  for  a fixed  false  alarm  l^te  (see  Anderson, 
1958).  Eecause  of  the  distributional  problems  associated  with 
the  output  of  this  quadratic  detection  filter  (QDF) , it  is  useful 
to  consider  several  versions  of  an  optimum  linear  filter.  We 
shall  consider  the  linear  detection  filter  (LDF)  which  gives  the 
highest  detection  probability  of  any  linear  filter  and  the  linear 
matched  filter  (LMF)  which  is  optimum  if  n is  a white  nois?  pro- 
cess with  equal  power  under  both  and  Matched  filters 

have  been  applied  to  detecting  long  period  signals  by  Glover  and 
Alexander,  1968,  and  Alexander  and  Lambert,  1971,  wno  used  a 
seismic  signal  as  the  reference  event.  Capon  et  al.,  1969,  pro- 
duced comparable  results  using  a chirp  waveform  instead  of  a 
seismic  signal. 

In  this  report  we  will  consider  the  application  of  the 
classical  spectral  measures  and  optimum  linear  and  quadratic 
filtering  to  the  population  of  earthquakes  and  explosions  shown 
in  Figures  Al  and  A2  of  the  Appendix.  For  referral  purposes  the 
discriminants  considered  are  listed  below: 

(1)  Linear  Detection  Filter  (LDF) 

(2)  Linear  Matched  Filter  (IMF) 
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(3)  Quadratic  Detection  Filter  (QDF) 

(4)  Integrated  Spectra  (.4-. 8 Hz) 

(4) '  Integrated  Spectra  (1.4-1. 8 Hz) 

(5)  Complexity 

(6)  Complexity  vs  Integrated  Spectra  (.4-. 8 Hz) 

The  distribution  theory  for  discriminants  (1),  (2),  (4)  and  (4) 
will  be  g'ven  which  allows  the  explosion  false  alarm  and  detec- 
tion probabilities  to  be  explicitly  calculated.  The  procedure 
for  using  the  above  discriminants  will  be  tested  by  using  a ran 
dom  subset  of  the  total  to  derive  the  relevant  parameters  and 
then  testing  the  derived  discriminants  on  the  remaining  events. 
Sections  2,  3 and  4 give  the  equations  for  the  various  proposed 
discriminants  and  summarize  the  results  of  applying  the  discri- 
minants to  the  full  suite  of  events.  Section  5 presents  the 
results  of  using  discriminants  derived  from  a randomly  drawn 
learning  set  tc  discriminate  on  the  remaining  presumed  "new" 
events.  In  Section  b we  discuss  the  physical  basis  for  the 
discriminants  as  a guide  to  the  expected  behavior  on  new  events 
and  give  suggestions  for  further  research. 


LINEAR  FILTERING 


In  this  section  we  will  investigate  linear  filtering  as  a 
technique  for  discriminating  between  hypotheses  and  as 
specified  in  equation  (1.1).  To  be  more  precise,  assume  that 
under  the  hypothesis  Hj , j = 1,  2,  the  T dimensional  observed 
vector  x'  = [x(0) , x(l),...,  x(T-l)]  is  a non-stationary  normal 
process  with  mean  signal  vector  s!  = [Sj(0),...,  Sj(T-l)]  and 
covariance  matrix 

Ej  (x-s ^.)  (x-Sj)  ’ = Rj  • C2-1) 

Since  the  noise  is  stationary  R.  is  assumed  to  contain  covariance 
functions  of  the  form 

R.(k-£)  - j1'2  e27ri(k_il)£P.(f)df  (2.2) 

J -1/2  J 

in  the  kta  row  and  £th  column.  Hence,  an  equivalent  statement 
of  the  hypothesis  Hj  is  that  x is  a normal  process  with  mean 
s.  and  noise  power  spectrum  Pj(f).  Define  the  difference  between 
tne  earthquake  and  explosion  vectors  as 

5 = Si  - s 2 • (2.3) 

We  seek  a linear  filter,  say  b'x,  to  use  for  classifying  x 

into  H i or  H7.  An  observation  is  classified  into  if  b'x  > c 

and  into  Hn  if  b'x  _<  c.  A simple  possibility  is  matching  the 

mean  difference  vector  5 with  the  data  vector  x producing  the 

matched  filter  output  5'x.  If  R-,  = R^  = o I-j.  with  1^  the  T x T 

identity  matrix,  the  probability  of  detecting  H 2 for  a fixed 

false  alarm  rate  will  be  maximized.  If  = R2  = *'e• 

p ( f ) = P7(f)  = P(f) , the  classical  linear  discriminant  function 
1 1 *- 

3 ' R x maximizes  the  detection  probability  as  shown  in  Anderson, 
1958. 
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If  Rj  f , Kullback,  1959,  hr.s  shown  that  for  a i ixed 
false  alarm  rate  a,  a discriminant  filter  of  the  form 


y(t1,t2)  = S'R^x 


maximizes  the  detection  probability  where 


R = *1*1  + *2*2 


(2.43 


(2.5) 


is  a weighted  combination  of  the  auto-correlation  matrices  and 

tl  and  t2  are  chosen  (subject  to  the  cmdition  that  R be  positive 
definite)  to  maximize  the  detection  probability. 


pd(°0  = pH  (accept  H- } = $ 


with 


d + $ * (a) a 


a2  = 6 ,R'1R.f"16 
J J 

d = 6,F16 


4>(x)  = / e’u  ^2  — . 

/2tT 


(2.6) 


(2.7) 


(2.8) 


(2.9) 


We  have  chosen  to  maximize  the  probability  of  detecting  an 
explosion  for  a given  false  alarm  rate.  A parallel  investiga- 
tion of  the  filters  which  maximized  the  earthquake  detection 
probability  for  a given  explosion  false  alarm  rate  gave  approxi- 
mately the  same  filters  and  performance  results.  Anderson  and 
Bahadur,  1962,  discuss  the  admissibility  of  estimators  having 
this  same  general  form.  The  general  procedure  for  using 
equations  (2.5) -(2.9)  is  first  to  choose  an  a and  then  to  search 
over  tj  and  t2  for  the  maximum  pd(a)  as  defined  in  equation  (2.6) 

The  difficulty  in  applying  equation  (2.6)  arises  from  the 
computational  effort  involved  in  the  matrix  calculations  (2.7) 
and  (2.8).  If  T x T (in  this  case  T = 256)  matrices  are 
involved,  even  the  Levinson  recursion  would  lead  to  excessive 


r 


* 


compute)  time  because  of  the  number  of  repetitions  necessary  to 
maximize  (2.6).  However,  a spectral  approximate  n method  based 
on  the  finite  Fourier  transform  (FFT)  can  be  employed  which 
yields  a great  saving  in  computer  time.  Spectral  approximations 
were  proposed  for  a problem  involving  only  spectra  by  Capon, 

1965,  and  we  re^er  to  Shumway  and  Unger,  1974,  for  the  conditions 
under  which  the  following  approximations  are  valid.  The  impor- 
tant condition  is  that  all  spectra  involved  in  the  equations  be 
strictly  positi”e. 

Suppose  that  the  FFT  of  the  function  g(t)  , t = 0,  1,...,T-1 
is  defined  as  usual  by 


where 


T-l 


(X)  - rin  l g(t)e-i*nt 

n n 


n = 0 


n 


(2.10) 

(2.11) 


with  f^  corresponding  to  f in  equation  (2.2)  and  X^  denoting 
the  angular  frequency.  Then,  approximations  to  (2.7)  and  (2.8) 
are 

. T-l  - P.fX  1 

(2.12) 


- T-l  2 P,(X) 

of.  I |D(xn)l2  r1— 2- 

J n=0  n P(Xn) 


and 


d ■ 1 


(2.13) 


n=0  P(Xn) 


where  D(Xr)  is  the  FFT  of  6(t),  Pj(XR)  is  the  spectrum  of  x 

f(»„)  = 'iW  * lzW'  t2-14) 


under  and 


Shumway  and  Unger,  1974,  show  that  T Oj  and  T Oj  have  the 
same  finite  limit  and  that  T_1d  and  T_1do  have  the  same  finite 
limit.  The  linear  detection  filter  (2.4)  (LDF)  may  then  be 
approximated  by 


i 


i 


i 

i 
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(2.15) 


T-l  D*(Xn)X(Xn) 
y Ct1,t2)  = I 


n=0  P(Xn) 

where  X(Xn)  is  the  FFT  of  the  data  trace  x(t)  and  * denotes  the 
complex  conjugate. 

We  refer  to  (2.15)  from  here  on  as  the  linear  detection 
filter  (LDF) . 


The  linear  matched  filter  ( LMF)  is  given  by 

Z = 5 'x  = j D*(Xn)X(Xn) 
r =0 

and  has  detection  probability  approximated  by 


Pa  (a)  = <t> 


d'  + $ *(01)0° 


(2.16) 


(2.17) 


with 


and 


(o|)2  - T |D(X  ) I^PjtX  ) 
J n = 0 J 


T-l 


d'  - Y |0(Xn)|2. 

n = 0 


(2.18) 


(2.19) 


In  order  to  apply  the  LDF  and  LMF  to  th ; full  suite  of  40 
earthquakes  and  26  explosions  given  in  the  Appendix,  estimates 
for  s^t),  s2(t),  PxUn)  and  P2(Xn)  are  needed.  I f we  regard 
x-k(t)  as  the  kth  observed  event  of  type  j,  j=l>2,  k=l,...,N^, 
we  may  estimate  the  j ^ signal  process  by 


N. 


(2.20) 


and  the  j*^  noise  spectrum  by 

N . 


Pj(x„)  - Ujk(xn)  - S,txn)i 


(2.21) 


-9- 


where  the  capital  letters  denote,  as  usual,  the  FFT.  These 
estimators  (after  smoothing  with  a three  point  running  average) 
are  shown  on  the  left  side  of  Figure  1.  In  this  example  T - 2 5 1 
points  are  sampled  at  10  points  per  second  yielding  a folding 
frequency  of  5 Hz. 


search  over  the  possible  values  for  tj  and  yielded 


1 = *2  ~ 1 as  the  solution  which  maximized  the  cetection  proha 


bility  for  a false  alarm  rate  for  explosions  of  a = .001.  The 
predicted  explosion  detection  probability  for  this  false  alarm 
rate  is  .997  with  a threshold  for  the  detector  established  as 
C = -8.51,  i.e.  we  accept  1^  (earthquake)  if  y°  (1,1)  > -8.51 
or  ll2  (explosion)  if  y°  (1,1)  < -8.51.  The  filter  for  this 
case  is 

T-l  D*  ( X ) X ( A ) 


y°CM)  - f FTrn~i~Mx  ) 

n=0  W VV 


so  that  we  may  look  at  it  as  a matching  of  the  signal  difference 
vectoi  LH  X^)  with  the  data  X(X^)  where  the  weights  are  inversely' 
proportional  to  the  average  spectrum.  The  vertical  axis  in 
Figure  2 shows  the  output  of  the  LDF  and  it  is  clear  that  38 
of  40  earthquakes  and  22  of  26  explosions  would  be  detected 
tor  a .001  explosion  false  alarm  rate.  Adjusting  the  threshold 
slightly  picks  up  one  more  explosion.  The  predicted  detection 
pi  oh  ah i lity  for  explosions  i>  .997.  Several  difficult  events 
are  causing  some  problems.  The  main  defect  appears  to  be 
the  assumption  that  each  earthquake  has  the  same  noise  spectrum 
and  that  each  explosion  has  the  same  noise  spectrum  even  though 
the  spectra  for  earthquakes  and  explosions  may  differ.  In  the 
discussion  we  trace  this  fact  to  the  possibility  that  some  of 
the  explosions  are  cratering  experiments,  and  that  some  of  the 
earthquakes  are  on  P-wave  radiation  nulls  which  are  more  effec- 
tive at  long  periods  than  at  short  periods.  The  missed  events 
are  marked  in  the  Appendix  with  a (1)  and  (2)  depending  on 
whether  the  LDF  or  LMF  failed  to  detect.  For  an  explosion  false 


TABLE  I 


Estimated  Filter  Parameters  for  the 
Full  Population  and  Sub  - Population* 


Filter  Full  Sub 


Characteristics 

LDF 

LMF 

LDF 

LMF 

Mean  (H^) 

1.0 

.01 

3.4 

.09 

Mean  (H-) 

-15.9 

-.59 

-13.5 

-.53 

Var  (HJ 

9 . 5 

.041 

9.8 

.042 

Var  f 

7.4 

.012 

7.0 

.011 

Threshold 

- 8.51 

-.32 

- 6.27 

- .39 

(1,1) 

(1,1) 

he  sub  - population  or  learning  set  contained  earthquakes 
H,  4,  6,  8,  10  , 11  , 13,  14  , 15  , 16,  21,  24,  25  , 26  , 27, 
.3,  30,  31,  32,  35,  36,  37,  40)  and  explosions  (2,  3,  6, 
7,  8,  9,  11,  15,  17,  18,  21,  23,  24,  25,  26). 


alarm  p*-obabilit'  of  .01  and  detection  probability  .99,  the 
linear  matched  .ter  (LMF)  shown  as  the  horizontal  axis  in 
Figure  2 detected  23  of  26  explosions  but  only  52  of  40  earth- 
quakes. The  filter  means  and  variances  are  shown  in  Table  I. 

We  note  here  that  the  proper  method  for  evaluating  the  detectors 
is  to  randomly  choose  a subset  of  events  to  estimate  the  filter 
parameters  and  then  apply  the  filters  derived  to  the  remaining 
events.  This  procedure  will  be  followed  in  the  Experimental 
Results  section. 


fl 
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OPTIMUM  QUADRATIC  FILTERING 


i 

i 


While  tiie  LDF  of  the  previous  section  is  the  optimum  linear 
filter,  it  is  not  tiie  optimum  filter  in  general  when  the  spectra 
of  the  earthquakes  and  explosions  are  different.  In  this  case, 
the  likelihood  approach  yields  a quadratic  form  for  determining 
whether  an  observation  x comes  from  or  U 2 (see  Anderson, 
1958).  In  this  case  the  log  likelihood  ratio 

S(x)  = -i(x-s1)  'R'Vx-sp  + -j(x-s2)  ,R21(x-s2)  (5.1) 


classifies  by  comparing  the  weighted  distance  between  x and  s 
with  the  weighted  distance  between  x and  s0.  A spectral  approxi- 
mation to  S(x)  is  given  by 


S°(x) 


1 V I^V'WI2  . 1 V |XCAn)-S2CA„)  I2 
- n=0  TO  ' * n*0  TO? 


(3.2) 


where  the  FFT's  anc  power  spectra  are  as  given  in  Section  2 and 
the  difference  T ^ |S(x) -S°(x) | converges  (wi'.h  probability  1) 
to  zero. 

One  difficulty  with  applying  S(x)  or  S°(x)  as  a discriminant 
is  that  the  probability  distributions  used  for  computing  detec- 
tion probabilities  are  not  tractable.  However,  it  is  known  that 
the  theoretical  detection  probability  of  the  quadiatic  detection 
filter  will  exceed  that  of  any  LDF.  The  only  real  problem  are 
setting  the  proper  threshold  and  worrying  about  the  effects  of 
the  estimated  spectra  and  signal  means  on  the  value  of  (5.2). 

Figure  5 shows  the  result  of  applying  the  discriminant  (3.2) 
to  the  full  suite  of  events.  We  see  that  a threshold  set  at 
about  . 10  would  have  detected  58  out  of  40  earthquakes  and  24 
out  of  2b  explosions.  The  misclassifications  are  marked  on  the 
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Figure  3.  Output  of  Quadratic  Detection  Filter  (QDF)  and  Linear 
Detection  Filter  (LDF)  Designed  and  Applied  to  Full  Suite  of 
I. vents.  Quadratic  threshold  at  0.1  estimated  visually. 


the  original  traces  in  Figures  A1  and  A2  with  a (3).  Again, 
while  the  i|ualitative  characterist  ics  of  the  QDF  appear  to  be 
excellent,  we  defer  making  any  comparative  statements  before 
the  1 xperinental  Results  section. 


SPECTRAL  RATIOS  AND  COMPLEXITY 


The  basic  measure  of  interest  here  is  the  integrated  spectrum 
as  measured  by  the  periodogram  accumulated  over  some  frequency 
banu  of  interest.  In  spectral  ratio  computations,  it  is  customary 
to  use  a low  frequency  band  as  the  range  .4-. 8 Hz  and  a high 
frequency  band  as  the  range  1.4-1. 8 Hz. 

Ihe  integrated  spectral  estimate  centered  at  f Hz  or  at 
radial  frequency  Aq  = 2 tt f q with  f a multiple  of  1/T  was  defined 
as  (L  < ( 1/2)T) 


LP(Xo) 


1/2(L-1) 
s=- 1/2 (L- 1) 


(4.1) 


where  X(.)  denotes  the  FFT  of  the  time  series  x(.).  If  the 
spectrum  is  relatively  smooth  o/cr  the  band  of  interest  then, 
approximately,  under 


LP(Ao) 


LP - ( A ) 
j v o' 


(4.2) 


\ I 7 

where  ~ denotes  "is  distributed  as"  and  X^L  denotes  a chi-square 
variable  with  2L  degrees  of  freedom.  Note  that  the  spectrum 
Pj(*0)  is  regarded  here  as  containing  both  a stochastic  signal 
and  a stochastic  noise  process.  Equation  (4.2)  allows  one  to 

A 

set  a threshold  for  the  sample  cumulative  spectrum  LP(Aq)  if  we 
know  the  earthquake  spectrum  Pj(A  ) and  the  explosion  spectrum 
P2^o^*  For  thc  cxamPle  suppose  that  P^Aq)  > P-,(>0)  and  we 
accept  H^ (earthquake)  if  LP(Aq)  > C and  we  accept  H 7 (explos ion) 
if  LP(i0)  i c-  Then,  the  explosion  false  alarm  rate  is 

° = ph1  |LF(a0)  - c(a)|  = p {#  I ri^(Ao)  > • (4-3) 


Solving  the  above  equation  for  C(a)  determines  the  threshold, 
with  the  detection  probability  given  by  (4.3).  If  there  arc 
spectral  estimates  in  another  non-overlapping  frequency  band 
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centered  at  A^,  the  classical  "spectral  ratio"  will  have  an  F 
distribution. 

Figure  4 shows  the  cumulative  spectra  for  two  band,  .4-. 8 Hz 

(f  = .6  Hz)  and  1.4-1. 8 Hz  (f  = 1.6  Hz).  Since  the  traces  are 
o o 

normalized  to  a unit  amplitude  which  tends  to  have  a period  near 
one-second,  the  low  and  high  frequency  spectra  tend  to  measure  the 
energy  content  relative  to  1 Hz.  From  the  plots  we  observe  that 
the  low  frequency  integrated  spectrum  discriminates  completely. 

As  we  shall  see  in  the  Discussion  section,  this  is  probably  due 
to  the  cancellation  of  P by  pP  at  low  frequencies  ^or  explosions. 
The  threshold  values  are  derived  by  noting  that  in  this  case 
I = 11  in  (4.2)  with  X^C  .01) /22  and  x|2  ( .99) /22  given  by  .434 
and  1.83  respectively.  In  order  to  arrive  at  reasonable  values 
for  LP1(Xq)  and  LP2(Xq),  the  integrated  spectra  were  averaged 
for  earthquakes  and  explosions  to  give  = .5802  and 

P0(X  ) = .1372  for  the  band  .4-. 8 Hz.  A cutoff  point,  set  at 
C(.01)  = .22  yields  a predicted  detection  of  Pd(.Ol)  = .99.  In 
this  case  the  sample  results  are  that  24  out  of  26  explosions 
and  32  out  of  40  earthquakes  are  detected.  Adjusting  the  thresh- 
old detects  3 additional  earthquakes.  The  high  frequency  band 
(1.4- 1.8  Hz)  performs  poorly  and  these  values  were  ignored.  The 
incorrect  classifications  for  the  spectra]  discriminator  are 
indicated  by  the  code  (4)  in  Figures  A1  and  A2. 

Complexity  is  measured  by  some  function  proportional  to  the 
total  integrated  power  in  the  signal  coda.  Here,  the  coda  is 
defined  as  the  part  of  the  signal  beginning  after  the  60th  point 
(6  seconds).  Figure  5 shows  complexity  plotted  against  the 
integrated  low  freauency  spectrum  and  we  note  that  with  a graphi- 
cally determined  cutoff  at  .25  we  detect  22  of  26  explosions  and 
36  out  of  40  earthquakes.  The  combining  of  complexity  and  inte- 
grated low  frequency  spectral  measurements  tends  to  improve  the 
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Figure  4.  Components  of  Spectral  Ratio  for  Full  Suite  of 


overall  discrimination  if  one  is  willing  to  draw  a negatively 
sloped  separation  line.  Such  a line  can  (see  Figure  5)  be 
adjusted  to  detect  23  of  26  explosions  and  3?'  of  40  earthquakes. 

A negatively  sloped  line  has  been  obtained  by  Anglin,  1971. 
Blandford,  1974,  has  also  discussed  the  simultaneous  increase 
in  complexity  and  decrease  in  low  frequency  power  for  earthquakes 
from  a theoretical  point  of  view. 

It  is  worth  noting  that  the  relative  ineffectiveness  of 
the  1.4-1. 8 Hz  band  energy  can  be  deduced  by  inspection  of 
Figures  15  and  14  in  LaCoss  (1969-24).  These  figures  plot 
0.4-0. 8 Hz  and  1.4-1. 8 Hz  energy  versus  m^ , and  one  sees  a far 
greater  separation  for  the  low  frequency  band. 
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EXPERIMENTAL  RESULTS 

All  of  the  discriminants  involved  here  need  some  prior  data 
in  older  to  establish  a decision  criterion  for  future  observations. 
1 he  filtering  discriminants,  LDF , LMF , QDF  all  need  a reference 
earthquake  and  a reference  explosion  to  match  against  a new  event 
of  unknown  origin.  In  addition,  the  LDF  and  QDF  need  prior  esti- 
mates of  the  earthquake  and  explosion  spectra  in  order  to  determine 
the  proper  weighting  functions  in  the  filter  and  to  establish  a 
threshold  value  for  the  discriminant  when  it  is  applied  to  a new 
data  trace.  The  "spectral  ratio"  (normalized  integrated  spectrum, 
.4-. 8 Hz)  and  complexity  also  require  an  initial  learning  set  to 
determine  a threshold  value  for  new  data. 

In  order  to  test  the  discriminants  under  realistic  conditions, 
a sub-population  of  23  earthquakes  and  15  explosions  was  drawn 
randomly  from  the  full  suite  of  40  earthquakes  and  26  explosions. 
The  mean  signals  and  spectra  for  the  subset  are  shown  on  the 
right  hand  half  of  Figure  1.  It  is  clear  that  *he  waveforms  and 
spectra  determined  for  the  random  sample  do  not  differ  greatly 
from  those  determined  from  the  full  suite  of  events.  Table  I 
compares  the  means,  variances  and  threshold  values  for  the  linear 
detection  and  linear  matched  filters.  All  characteristics  notice- 
able in  the  full  set  seemed  to  be  present  in  the  learning  set. 

The  parameters  given  by  the  learning  set  were  then  used  to 
derive  the  discriminants  to  be  applied  to  the  events  not  in  the 
learning  set,  now  presumed  to  be  new  events  of  unknown  origin. 

Table  II  summarizes  the  results  of  applying  the  discriminants 
derived  from  the  learning  set  to  both  the  learning  population 
and  the  population  of  new  observations.  The  cutoff  values  deter- 
mined from  the  sub -population  for  the  quadratic  detector  spectral 
ratio  and  complexity  were  .10,  .22  and  .25  respectively.  Several 
qualitative  observations  can  be  made  about  the  various  detectors, 
lhe  quadratic  filter  had  a perfect  record  on  the  learning  set 


TABLE  II 

Summary  of  Experiment  Using  the 
Linear  Detection  Filter  (LDF)  , Linear  Matched  Filter  (LMF) , 
Quadratic  Detection  Filter  (QDF) , Spectral  Ratio  (SR)  , 
Complexity  (C) , and  Complexity  vs  Spectral  Ratio  (CSR) 


Learning  S *t 


Total 

Number 

Detected 

LDF 

LMF 

QDF 

SR 

C 

CSR 

Number  of 
Earthquakes 

* 23 

22 

23 

23 

20 

21 

23 

Number  of 
Explosions 

15 

14 

11 

15 

14 

12 

13 

Tes  t 

Set 

Total 

Number 

Detected 

LDF 

LMF 

QDF 

SR 

C 

CSR 

Number  of 
Earthquakes 

17 

17 

15 

15 

13 

15 

15 

Number  of 
Explos ions 

11 

10 

9 

9 

10 

10 

10 
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while  making  a total  of  4 misclassif ications  in  the  test  set. 

This  might  indicate  a sensitivity  to  the  assumption  that  the 
residual  spectra  for  earthquakes  taken  together  and  explosions 
taken  together  are  identical.  The  spectral  ratio  and  complexity 
performed  about  equally,  with  the  complexity  measure  doing 
slightly  better  on  the  test  set.  The  linear  detection  filter 
(LDF)  had  the  best  overall  performance  record  making  2 misclassi- 
fications  in  the  learning  set  and  only  1 misclassification  in  the 
test  set.  The  linear  matched  filter  was  roughly  comparable  to 
the  quadratic  detection  filter,  spectral  ratio  and  complexity, 
each  of  which  had  3 to  5 misclassifications  in  the  test  set. 

The  above  results  are  based  on  a relatively  small  learning 
set  and  a relatively  small  test  set.  Hence,  the  slight  superio- 
rity of  optimum  filtering  is  not  of  great  importance.  The  assump- 
tions which  guarantee  that  optimum  filtering  be  optimum  are  only 
approximately  satisfied.  Variations  are  present  in  the  individual 
error  spectra  and  the  use  of  estimated  parameters  may  cause  per- 
formance fluctuations  in  new  data  samples.  However,  the  results 
are  encouraging  enough  to  speculate  that  optimum  filtering  may 
be  as  effective  as  the  classical  discriminants  when  applied  to 
large  data  bases. 


-24- 


i 


i 


r 


DISCUSSION 


The  explosion  which  is  most  earthquake- like  is  number  19. 

This  event  is  in  the  region  of  the  proposed  Pcchora-Kama  Canal 
(Nordyke , 1975)  where  several  excavation  experiments  have  been 
proposed.  It  seems  quite  likely,  therefore,  that  this  event 
was  an  excavation  shot.  Thus  no  pP  would  have  been  generated, 
the  low  frequency  energy  would  not  be  canceled,  and  the  overall 
low- frequency  appearance  of  the  record  would  be  explained. 

Similarly,  Marshall  (1972)  has  reported  that  cratering 
experiments  have  been  carried  out  near  the  Kazakh  test  site. 

The  similarity  of  explosion  18  to  explosion  19  suggests  that  it 
also  is  a cratering  experiment,  thus  also  explaining  its  failure 
to  discriminate.  Discrimination  failures  for  cratering  ex plo- 
sions arc  not,  of  course,  a serious  matter  in  the  context  of  a 
test-band  treaty.  If  these  two  explosions  are  excluded,  one 
may  say  that  there  are  no  explosion  discrimination  failures  for 
the  0.4-0. 8 spectral  ratio. 

One  might  have  expected  a failure  for  explosion  22,  reported 
by  Marshall  (1972)  to  have  been  at  an  announced  depth  of  1.53  km, 
and  with  an  apparent  pP-P  time  of  1.2  seconds.  With  such  a large 
delay  the  0.4-0. 5 Uz  energy  would  not  be  severely  reduce'  by  the 
spectral  null  at  zero.  However,  the  null  at  0.8  Hz  apparently 
eliminates  enough  low-frequency  energy  to  enable  the  0.4-0. 8 Hz 
spectral  discriminant  to  work.  Thus  the  explosion  is  discrimi- 
nated despite  the  ineffectiveness  of  the  null  at  zero  which  is 
the  fundamental  physical  phenomenon  underlying  the  discrimina- 
tion of  the  other  events. 

Another  explosion  was  set  off  May  21,  1908  at  Bukara 
(Marshall,  1972)  at  an  announced  depth  of  2.45  km.  Marshall 
obtains  a delay  oi  l.i>2  seconds  which  would  place  spectral  nulls 
at  0.62  and  1.24  Hz.  We  find  that  this  event  also  classifies  as 
an  explosion  on  the  0.4-0. 8 Hz  spectral  ratio  discriminant.  It 
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would  therefore  be  difficult  to  evade  this  discriminant  by  setting 
off  the  component  explosions  of  a shot  array  at  considerable  depth. 
If  success  were  possible  at  some  depth,  one  would  then  have  to 
use  negative  discrimination  on  long-period  (10-20  sec)  P. 

We  must  keep  in  mind  that  there  is  every  reason  to  expect 
an  occasional  earthquake  which  will  classify  as  an  explosion  on 
the  0.4-0. 8 Hz  spectral  discriminan t . A shallow  (1-2  km),  small 
45°  dip-slip  earthquake  would  send  out  identical  P-waves  up  and 
down,  resulting  in  P-pP  cancellation.  These  are  the  same  events 
which  can  show  compres sional  first  motions  at  all  teleseismic 
distances;  and  which  are  difficult  to  discriminate  by  Ms:mb, 

Douglas  et  al.  (1973),  Blandford  (1974).  Another  possibility  is 
that  we  may  find  ourselves  on  a node  of  the  1 radiation  pattern, 
and  that  this  pattern  will  be  less  smoothed  out  by  inhomogeneities 
at  long  periods  than  at  short  periods. 

In  a parallel  study  we  have  applied  a more  severe  test  of 
the  idea  that  low  frequencies  are  lacking  in  underground  explo- 
sions by  applying  an  0.5-0. 5 Hz  passband  filter  to  the  events  in 
Appendix  I.  Several  events,  including  3 of  the  5 earthquakes 
which  fail  the  0.4-0. 8 Hz  spectral  discriminant,  show  no  detect- 
able signal  in  the  time  domain.  One  can  say  that  the  minimum 
level  is  greater  than  that  of  any  explosion.  However,  two  of  the 
earthquakes,  numbers  22  and  39,  still  classify  as  explosions. 

The  nature  of  these  events  is  unknown  at  the  present  time.  If 
we  may  assume  that  we  have  a fair  sample  of  earthquakes,  then 
at  single  stations  we  must  expect  a false  alarm  rate  of  5-151 
on  new  earthquake  data.  If  a substantial  portion  of  these  false 
alarms  are  due  to  the  radiation  pattern  problem  discussed  above, 
then  their  number  would  be  substantially  reduced  if  observations 
were  available  from  a few  well-distributed  stations.  A suitable 
method  of  analysis  might  be  to  define,  in  the  conventional  way, 

an  m,  measured  on  traces  filtered  through  a low-frequency  pass- 

fa  . 

band,  and  to  discriminate  on  the  difference  between  this  m^  and 

the  conventional  one. 


-26- 


— - 


Together  with  these  results  we  have  "0"$  false  alarm  rate 
for  explosions  at  "reasonable"  and  even  large  depths  and  a 
100°o  false  alarm  rate  for  cratering  explosions. 

This  discriminant  would  be  extremely  difficult  to  spoof  with 
a shot  array.  Each  individual  shot  would  have  less  low-frequency 
energy  than  required.  As  we  go  to  lower  frequencies  where  the 
low- frequencies  can  be  superposed  while  the  shorter  frequencies 
cancel;  the  amount  of  low-frequency  energy  steadily  decreases. 
Overburying  could  possibly  be  a solution  if  care  is  taken  to 
account  for  the  behavior  of  all  the  spectral  nulls,  see  the 
discussion  above  of  explosion  22.  It  seems  conceivable  that  a 
properly  timed  array  could  sufficiently  cancel  the  predominant 
energy  in  the  waveform  while  not  cancelling  the  low-frequency 
energy,  but  the  design  analysis  would  have  to  be  much  more 
subtle  than  in  the  corresponding  M^m^  ose,  and  long-period 
(10-20  sec)  P of  possibly  detectable  for  an  earthquake  or  com- 
parable sice  would  still  remain  as  a worry  to  the  evader  from 
the  point  of  view  of  negative  discrimination. 

Future  research  should  concentrate  on  the  details  of  measure- 
ment of  the  low  frequency  P-wave,  long-period  P discrimination 
and  thresholds,  applicability  of  the  low-frequency  P discriminant 
as  a function  of  distance,  compilation  of  an  extensive  data  base 
of  earthquakes  to  "regionalize"  the  discriminant,  and  investiga- 
tion of  the  depth  and  fault-plane  solutions  for  those  earthquakes 
which  classify  as  explosions. 
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APPENDIX  I 

Basic  Data  Traces  with  Misclass if ied  Events 
in  tiie  Full  Population  Discriminants  Noted  as: 

(1)  Misclassi f ied  by  linear  detection  filter  (LDF) 

(2)  Misclassi fied  by  linear  matched  filter  f LMF) 

(3)  Misclassi fied  by  quadratic  detection  filter  (QDF) 

(4)  Misclassified  by  spectral  ratio  (Int.  .4-. 8 Hz) 

1.5)  Misclassified  by  complexity 

(.0)  Misclassified  by  complexity  vs  low  frequency  spectrum 


Seismogram  numbers  arc  given  beneath  the  event  numbers. 
Numbers  in  parentheses  after  certain  event  numbers  are  for 
reference  to  earlier  studies  using  tnese  events.  The  source 
for  event  parameters  arc  given  with  each  data  trace.  All 
events  near  50°N,  78°E  are  presumed  explosions.  The  other 
presumed  explosions  have  been  identified  as  such  by 
Nordvke  (.1975)  or  Marshall  (1972)  citing  USAEC  or  SIPRI 
sources.  the  order  of  event  parameters  is  date,  latitude, 
longitude,  , and  depth. 
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APPENDIX  II 

Writeups  of  Relevent  Computer  Programs 


■ 
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- 1- 


A.  IDENTIFICATION 

1.  Title:  SPECES 

2.  Focus  Identification: 

3.  Category: 

4.  Entry  Points:  SPECES 

5.  Software  Systems:  Fortran  IV 

6.  Computer  Configuration:  IBM  S/360 

7.  Programmer:  Robert  Shumway 

8.  Contributing  Organization:  Geotech 

9.  Program  Materials  Available:  Source  Deck,  Listing 

10.  Date:  May  1974 


B.  PURPOSE 

Estimates  the  mean  signal  waveform  and  noise  spectrum  for  each 
of  a number  of  groups  of  time  series.  Each  series  is  adjusted  to 
unit  maximum  peak  to  peak  amplitude  and  zero  mean.  Estimates  for  the 
integrated  spectra  of  each  trace  in  two  bands  (f^,  ii)  and(f3,  f4) 
are  computed.  An  estimate  for  the  noise  spectrum  of  each  group  is 
calculated  by  subtracting  the  mean  signal  for  that  group  from  each 
trace  and  averaging  the  spectra  of  the  residuals. 

C.  USAGE 

1.  This  is  a Fortran  IV  main  program  which  calls  the  sub- 
routines DOT,  ADD,  CR,  FFT , PLOP,  RC , MINMAX,  SPLOT, 

REMAV , CUM,  RAVG , in  the  deck  and  ERASE,  CLOCK 

2.  Parameters 


Col. 

Parameter 

Descript  ion 

Card  1 

HIT 

Card  2 

1 

* 

Plot  character 

1-S 

NPT 

No.  of  pts.  in  each  extended  (2 

(1615) 

6-10 

II 

series 

No.  of  groups 

11-15 

NDATA 

No.  of  pts.  in  series  as  read 

16-20 

ISM 

No.  of  pts.  over  which  spectrum 

Card  3 

1-10 

SR 

smoothed 

Sampling  rate  in  pts/sec. 

TTFT072) 

11-20 

R 

Gain  adjusted  to  amplitude  R 

Card  4 

1-10 

FR1 

First  frequency  interval 

(4P10.2) 

11-20 

FR2 

Cumulated  is  (FR1,  FR2)  Hz 

21-30 

FR3 

Second  frequency  interval 

31-40 

FR4 

Cumulated  is  (FR3,  FR4)  Hz 

An- a 


Card  5 
1015 


Col.  Parameter 


1-5 

6-10 


N (1 ) 


Description 

Number  of  time  series  in  group  I 
is  N ( I ) 


Card  6 (Etc)  1-5 


NCI  I) 

NSKP 


One  less  than  desired  start  point 
for  each  series 


Space  Required:  72K 

Alarms : None 

Error  Returns : None 

Tape  Mountings:  Input,  SYS005,  Plot,  SYS004 

Formats:  See  parameters 

External  Symbols  and  Subroutines: 

Selective  Jump  and  Stop  Settings:  None 

Timing:  13  min. 

Cautions : None 

References:  SDAC  Report,  An  Examination  of 

Some  New  and  Classical  Short  Period 
Discriminants . 
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A.  IDENTIFICATION 

1.  Title:  LINDIJ 

2.  Focus  Identification: 

3.  Category: 

4.  Entry  Points:  LINDI1 

5.  Software  Systems:  Fortran  IV 

6.  Computer  Configuration:  IBM  S/360 

7.  Programmer:  Robert  Shumway 

8.  Contributing  Organization:  Geotech 

9.  Program  Materials  Available:  Source  Deck,  Listing 

10.  Date:  May  1974 


B.  PURPOSE 

Searches  the  two  group  admissible  linear  filters  for  the  one 
that  maximizes  the  detection  probability  as  a function  of  several 
input  false  alarm  rates.  The  mean  signals  and  noise  rpectra  Pj  and 
?Z  for  each  of  two  groups  (say,  from  SPECES)  are  required  as  input. 
The  means  and  variances  and  threshold  values  for  the  linear  detec- 
tion and  linear  matched  filters  are  computed.  The  detection  pro- 
babilities are  plotted  as  a function  of  the  weight  parameters 
applied  to  the  two  group  spectra.  Two  weighting  schemes  P1-0P2 
and  P2-0P1  are  considered  as  a function  of  0. 


C.  USAGE 

1.  This  is  a Fortran  IV  main  program  which  calls  the  sub- 
routines CR,  FFT , PLOP,  RC,  MAD,  MINMAX,  PROBF , QUANTF , 
SPLOT  ii  the  program  as  well  as  CLOCK,  ERASE,  SQRT,  LOG, 
PLOT. 

2.  Parameters 


Col. 

Parameters 

Description 

Card  1 

1 

* 

Plot  Character 

tait 

Card  2 

1-5 

NPT 

No.  of  pts.  (2^) 

(1015) 

6-10 

NIT 

No.  of  iterations 

11-15 

MIT 

Method  of  iteration, if  1 over  the 

16-20 

NALP 

admissible  thetas,  if  2 over  thetas 
read  in 

No.  of  false  alarm  rates  considered 

A-H-1 
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Col . 

Parameter 

Description 

Card  3 
(6F10  2) 

1-10 

R\NGE 

Range  of  theta  searched  if  MIT  is 

Card  1 

1-10 

THE (1 ) 

THE(I)  is  the  I^h  value  of  theta 

(6  FI  0 . 2 ) 

THE (NIT) 

if  MIT  is  2.  If  MIT  is  1 , leave 
those  cards  out. 

Card  5 
(4F20. 5) 

1-20 

ALPHA  (1) 
ALPHA (NALP) 

False  alarm  rate  I is  ALPHA(I) 

Card  6 (Etc) 
(8F10.3) 

1-80 

FI  (1) 

FI (NPT) 
F2  (1) 

F2 (NPT) 

Mean  Signal  of  1st  Group 
Mean  Signal  of  2nd  Group 

Card  7 (Etc) 
(8E10.3) 

1-80 

FI  (1 ) 

• 

FI (NPT) 

F2  (1 ) 

• 

F2 (NPT) 

Noise  Spectrum  of  1st  Group 
Noise  Spectrum  of  2nd  Group 

3.  Space  Required:  58K 

4 . Alarms : 

None 

5.  Error  Returns:  None 

6.  Tape  Mountings:  Plot,  SYS004 

7.  Formats:  See  parameters 

8.  External  Symbols  and  Subroutines: 

9.  Selective  Jump  and  Stop  Settings:  None 

10 . T iming : 30  min . 

11.  Cautions:  None 

12.  References:  SDAC  Report,  An  Examination  of 

Some  New  and  Classical  Short  Period 
Discriminants . 
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A.  IDENTIFICATION 


1.  Title:  LINDI2 

2.  Focus  Identification: 

3.  Category: 

4.  Entry  Points:  LINDI2 

5.  Software  Systems:  Fortran  IV 

6.  Computer  Configuration:  IBM  S/360 

7.  Programmer:  Robert  Shumway 

8.  Contributing  Organization:  Geotech 

9.  Program  Materials  Available:  Source 

10.  Date:  May  1974 


Deck,  Listing 


B.  PURPOSE 

Applies  the  linear  matched  filter  and  a collection  of  linear 
detection  filters  indexed  by  the  spectral  weighting  parameter 
0 ( L P 1 ~ ©P 2 ] or  [P2’0Pl])  to  a sample  of  time  series.  The  new 
series  are  classified  into  two  groups  using  the  input  critical 
values . 


C.  USAGE 

1.  This  is  a Fortran  IV  main  program  which  calls  the  subroutines 

T1  PT  tm  r\n  /\m  _____ 


ERASE,  PLOT 

• 

2. 

Parameters 

Col. 

Parameter 

Card  1 

1 

* 

(Al) 
Card  2 

1-5 

LF 

(1015) 

6-10 

NTHETA 

11-15 

NE 

16-20 

NDATA 

Card  3 

1-20 

THE(l) 

(4F20.5) 

21-40 

Card  4 

THE (NTHETA) 

1-10 

ID(1) 

(1015) 

11-20 

Description 
Plot  Character 

No.  of  pts.  in  extended  series  (2N) 
No.  of  values  of  0 to  be  considered 
No.  of  time  series  to  be  classified 
No.  of  data  points  read  in 
Values  of  parameter  0 


ID (NTHETA) 


Identification 
desired  and 


is 


which  is  1 if  Pj-OP-. 

2 if  P2*0Pi  is  desired 


A -I  i-  *• 


L. 


t*,  j"  m 


1 

1 

- 6 

Col. 

Paramete  r 

Description 

Card  5 
'(8  FI  0.3) 

1 

Card  6 (Etc) 
T8T10.3) 

1-10 

11-20 

1-10 

11-20 

C(l) 

C(NTHETA) 

XM1(1) 

• 

Input  critical  values  corresponding 
to  the  9's.  The  last  O is  a dummy 
value  corresponding  to  the  linear 
matched  f i Iter  . 

Mean  Signal  for  group  1 

XMl(LF) 
XM2  ( 1 ) 

Mean  Signal  for  group  2 

Card  7 (Etc) 
(8E10.3) 

1-10 

11-20 

XM2CLF) 
FI  (1) 

Noise  spectrum  for  group  1 

Fl(LF) 
F2  (1) 

Noise  spectrum  for  group  2 

Card  8 (Etc) 
TTST~ 

1-5 

• 

F2(LF) 

NSKP 

Unc  less  than  the  start  point  for 
the  series  read  in 

3.  Space  Required:  57K 

4.  Alarms:  None 

5.  Error  Returns:  None 

6.  Tape  Mountings:  Input,  SYS005,  PK  YS004 

7.  Formats:  See  parameters 

8.  Fxternal  Symbols  and  Subroutines: 

9.  Selective  Jump  and  Stop  Settings:  None 

1 0 . Timing:  15  min, 

1 1 . Caut i ons : None 

12.  References:  SDL  Report,  An  Examination  of 

Some  New  and  Classical  Short  Period 
Discriminants . 
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A.  IDENTIFICATION 


1.  Title:  QUAD1 

2.  Focus  Identification: 

3.  Category: 

4.  Entry  Points:  QUAD1 

5.  Software  Systems:  Fortran  IV 

6.  Computer  Configuration:  IBM  S/3b0 

7.  Programmer:  Robert  Shumway 

8.  Contributing  Organization:  Geotech 

9.  Program  Materials  Available:  Source  Deck,  Listing 

10.  Date:  May  1974 


B.  PURPOSE 

Computes  a number  of  non-linear  discriminants  for  a collection 
of  time  series  including  (1)  Quadratic  detection  filter  output 
(2)  Spectra  accumulated  over  the  bands  (flf  f2)  and  (f3,  f ,)  , (3) 
Complexity,  (4)  Skewness  (third  moment),  (5)  kurtosis  (fourth  moment) 


C .  USAGE 


1.  This  is  a Fortran  IV  main  program  which  calls  the  sub- 
routines CR,  FFT , PLOP,  RC,  ADD,  SPLOT,  M INMAX,  REMAV , 
DOT,  CUM,  DIST,  RAVG , SK. 


2.  Parameters 

Col.  Parameter 


Card  1 

1 

* 

TaTT" 

Card  2 

1-10 

EP 

( 8 F 1 0 . 3 ) 

11-20 

FR1 

21-30 

FR2 

31-40 

FR3 

41-50 

FR4 

Card  3 

1-5 

LF 

1015) 

6-10 

NDATA 

11-15 

NE 

16-20 

1ST 

Description 
Plot  Character 

Small  value  to  add  to  each  noise 
spectrum 

Low  frequency  integrated 
Spectrum  (FR1,  FR2 ) 

High  frequency  integrated 
Spectrum  (FR3,  FR4) 

Length  of  extended  data  series  (2N) 
No.  of  data  pts.  read  in 
Number  of  time  series 
Complexity  is  computed  starting  at 
pt.  1ST  + 1. 
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Col . 


Parameter 


Pesc  r ipt  ion 


Card  4 (Etc)  1-10 
(8F10.3)  : 


Card  5 (Etc)  1-10 
(8I'10.3) 


Card  b (Etc)  1-5 

TTsl 


a 

4 

5 

6 

7 

8 
9 

10 

11 

12 


XM 1(1) 

Mean  signal  for  1st  type  of 

event 

XM1 (LF) 

XM2 (1 ) 

Mean  signal  for  2nd  type  of 

event 

XM2 (LF) 

fid) 

Noise  spectrum  for  1st  type 

of 

FI (LF) 

event 

F2  ( 1 ) 

Noise  spectrum  for  2nd  type 

of 

F2 ( LF) 

event 

NSKP 

One  less  than  the  start  point  for 

.:  64K 

each  series  read  in 

None 

: Input , 

SYS005,  Plot,  SYS004 

Alarms:  None 

Error  Returns : 


Formats:  See  parameters 

External  Symbols  and  Subroutines: 

Selective  Jump  and  Stop  Settings:  None 

Timing:  23  min . 

Cautions:  None 

References:  SDAC  Report,  An  Examination  of 

Some  New  and  Classical  Short  Period 
Discriminants . 
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